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The measurement of WZ production is an important 
test of the standard model (SM) of particle physics. WZ 
pairs are produced in both s-channel (qq — > W* —>■ WZ) 
and in t-channel (qq —> WZ) interactions. The WZ pro- 
duction is unique in that the s-channel mode of produc- 
tion provides sensitivity to the WW Z vertex which is 
governed by triple-gauge boson couplings (TGCs); the 
presence of anomalous couplings [1] could be an indi- 
cation of new physics at a higher mass scale leading to 
different rates and kinematic distributions than predicted 
by the SM. Furthermore, this process is an essential back- 
ground for Higgs boson searches at particle colliders be- 
cause the WZ decay into leptons is the primary back- 
ground to high mass Higgs boson searches in the three- 
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lepton signature, as well as an important background pro- 
cess in two- lepton Higgs boson analyses [2j. 

This Letter reports a measurement of the WZ pro- 
duction cross section and limits on anomalous TGCs us- 
ing a final state consisting of three charged leptons and 
one neutrino in pp collision data collected by the CDF 
II detector from 7.1 fb" 1 of integrated luminosity. The 
WZ — » Ivll decay, where I is an electron or muon and v is 
a neutrino contributing missing energy, allows for the re- 
construction of a variety of kinematic quantities that are 
utilized to distinguish signal from background by train- 
ing a neural network The cross section is then ex- 
tracted from the output of the neural network by fitting 
the shape with a maximum likelihood method. Limits 
on anomalous TGCs are set by analyzing the shape of 
the component of the momentum of the Z boson that is 
transverse to the beamline as a discriminant. 

The SU{2) L x U(1) Y part of the SM implies that W 
and Z weak vector bosons may interact via trilinear or 
quartic vertices. The WZ cross section is proportional to 
the interaction coupling strength predicted by the SM as 
e cot Ow i where e is the electric charge and 9w is the weak 
mixing angle. We want to measure the cross section of 
the process pp — > WZ, whose expected value in the limit 
of zero W and Z boson decay widths is 3.50 pb [!, 
We therefore measure the cross section by performing a 
likelihood fit to a parameter representing the ratio of the 
measured to expected WZ cross section, to be discussed 
below in this Letter. 

The cross section of this process was first reported by 
CDF with 1.1 fb- 1 1 6] in 2007. Subsequently, D0 re- 
ported with 1.0 fb" 1 0| in 2005, updated with 4.1 fb" 1 [8] 
in 2007, and again with 8.6 fb" 1 9] in 2012. Limits on 
anomalous TGCs were reported previously by D0 (loj . 

In the CDF II detector [llj, a particle's direction is 
characterized by the azimuthal angle <f> and the pseu- 
dorapidity 77 = — ln[tan(0/2)], where 9 is the polar an- 
gle measured with respect to the proton beam direction. 
The transverse energy Et is defined as Esin9, where E 
is the energy in the electromagnetic (EM) and hadronic 
calorimeter towers associated with a cluster of energy de- 
position. The transverse momentum pt is the particle's 
momentum component transverse to the beam line. The 
magnitude of the px for an electron is scaled according 
to the energy measured in the calorimeter in order to ac- 
count for momentum loss from final state radiation and 
bremsstrahlung. The missing transverse energy vector 
iftrp is defined as —'}2, i E T n T , where the index i loops 
over all calorimetric towers and n T is the unit vector in 
the transverse plane pointing from the interaction point 
to the energy deposition in calorimeter tower i. The $ T 
is corrected for the px of muons, which do not deposit 
all of their energy in the calorimeter, and tracks which 
point to uninstrumented regions in the calorimeter. The 
scalar missing transverse energy is defined as \$ T \ and 
denoted as $ T . Strongly interacting partons produced 



in the pp collision undergo fragmentation that results in 
highly collimated jets of hadronic particles. Jet candi- 
dates are reconstructed using the calorimeter signals and 
are required to have Et > 15 GeV and \rj\ < 2.5. Iso- 
lated lepton candidates are accepted out to an of 2.0 
for electron candidates and \rj\ of 1.0 for muon candidates. 

The experimental signature for the decay WZ — > Ivll 
is reconstructed as three charged leptons (electrons or 
muons) and $ T from the neutrino(s) that escaped unde- 
tected. Events are also detected if the W or Z decays to 
tau lepton(s) and those tau(s) subsequently decay to de- 
tectable electrons or muons-these events are considered 
part of the signal. Consequently, events containing three 
charged leptons, not all with the same charge, are se- 
lected from the data sample. The online event triggering 
and selection of lepton candidates are identical to those 
used in the search for SM Higgs bosons decaying to two 
W bosons at CDF Our baseline event selection is 
to require the leading-lepton's Et (or pr for muons) to 
be above 20 GeV (GeV/c) to satisfy the trigger require- 
ments, while the second and third leptons are allowed to 
have an Et (pt) as low as 10 GeV (GeV/c). Additionally, 
because the neutrino in the W — > Iv decay carries unde- 
tected energy, the WZ — > Ivll process tends to produce 
events with higher missing energy than the background 
processes-aside from it whose $ T distribution is also sim- 
ilarly high-valued, but is a nearly negligible background. 
We therefore require $ T > 25 GeV. 

Lastly, the dominant background remaining after these 
cuts is SM ZZ production. This motivates two more cuts 
to require one and only one Z 11 candidate in the event. 
We make a standard Z boson identification cut by requir- 
ing events to have a pair of same-flavor, opposite-signed 
leptons whose two-lepton mass falls within a window of 
±15 GeV/c 2 around the Z mass. This removes most of 
the SM backgrounds with no Z in the final state. We 
note that this cut reduces Zj — > llj events as the dilep- 
ton mass would not reconstruct back to the Z when the 
7 is emitted by one of the two leptons. 

To further reduce the ZZ — > llll background, we re- 
ject any event with an extra track with pt > 8 GeV/c, 
thereby rejecting events that may have a fourth lepton 
that failed to be identified. This cut reduces the remain- 
ing ZZ background by ~ 36% while leaving the WZ 
signal contribution essentially unchanged. Even so, ZZ 
remains the primary background in this measurement. 

There are several SM processes that result in a similar 
final state to WZ and are backgrounds in this measure- 
ment. The aforementioned ZZ — > llll process appears 
as a background when one of the four leptons fails to be 
reconstructed by the detector. This leaves three recon- 
structed leptons with the one lepton failing reconstruc- 
tion providing the missing energy signature. Drell-Yan 
events produced in association with hadronic jets that 
mimic the signature of a third lepton as well as Drell-Yan 
pairs produced with an associated photon that converts 
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TABLE I: Expected number of signal (WZ) and background 
events along with the total number of expected and observed 
events in the data. Uncertainties include all systematic un- 
certainties described in the text. 



Process 


Events 


ZZ 


3.6 


± 


0.5 


Z+jets 


3.4 


± 


0.8 


Z 1 


0.8 


± 


0.3 


tt 


0.1 


± 


0.04 


Total background 


8.0 


± 


1.0 


WZ 


47.4 


± 


4.8 


Total expected 


55.4 


± 


4.9 


Data 




64 





to an electron-positron pair via interaction with the de- 
tector are also significant backgrounds. Lastly, top quark 
pair production (tt — > W + bW~b) provides a minor con- 
tribution to the background when one of the subsequent 
&-quark jets mimics a lepton signature. The sum of these 
four backgrounds is quite small compared to the expected 
signal in the signal kinematics region. 

Generation of the background modeling-with the ex- 
ception of the Z+jets background-are determined by 
simulation. Events from WZ, ZZ, and tt are simu- 
lated using the PYTHIA [3] generator. The Z7 back- 
ground normalization is determined using the genera- 
tor described in Ref. [T3 |. The response of the CDF 
II detector is modeled with a GEANT3-based simula- 
tion program. The expected yields for each pro- 
cess are normalized to the cross sections calculated at 
partial next-to next-to-leading order (tt next-to- 
leading order (WZ and ZZ |4j]), or leading-order with an 
estimated normalization correction to account for higher- 
orders (Z7 [H|])- Efficiency corrections for the simulated 
detector response to lepton candidates are determined 
using samples of observed Z — > events. The Z+jets 
background normalization is calculated using the prob- 
ability that a hadronic jet will be reconstructed as a 
lepton candidate (the same as what is done in CDF's 
H — s- WW search [2]), which is measured in indepen- 
dent jet-triggered data samples. These probabilities are 
applied to the jet in the Z+jets data sample to estimate 
the number of such events that will pass the lepton identi- 
fication and signal selection criteria. The expected signal 
and background contributions are given in Table U along 
with the observed number of events. 

The dominant systematic uncertainties on the esti- 
mated contributions come from the luminosity measure- 
ment (6%) [l7| and the simulated acceptances of the sig- 
nal and background processes. The acceptance uncer- 
tainty due to the parton distribution function modeling 
ranges from 2.1% to 2.7% for the various processes. A 
10% uncertainty is assigned to WZ and ZZ processes 
for the kinematic differences between leading-order and 
higher-order calculations. The cross section uncertainty 



is 6% on the ZZ process, 7% on ti, and 5% on Z7. The 
Z"f process has another 20% uncertainty that accounts 
for possible mismodeling of the rate at which the 7 is 
misidentified as a lepton. Similarly, there is a 24.5% 
(23%) uncertainty on Z+jets (tt) for mismodeling the 
rate at which light jets (&-jets) are misidentified as a lep- 
ton. The uncertainty on the modeling of lepton identi- 
fication is 2% and of trigger efficiencies is 5.4%. Lastly, 
uncertainties on overall rates for the modeling of jets ac- 
counts for 1.2%. 

Within the signal kinematic region, we seek to further 
isolate the signal from background by utilizing a Neu- 
roBayes neural network treatment [3|. In general, the 
benefits of using a neural network (NN) over a simple 
counting experiment are two-fold: it can better isolate 
the signal from the background and provide a single dis- 
tribution from which the cross section value can be ex- 
tracted by fitting the data to the shape of the expected 
physical processes. We train a neural network with a 
combination of background events and simulated signal 
events. The input variables for the NN are kinematic 
quantities selected to exploit differences between signal 
and background distributions. Starting with many quan- 
tities that show relatively small differences in the dis- 
tributions of backgrounds and signal, a neural network 
will assign a numerical score whose distribution for back- 
grounds and signal will be better separated than in any 
single input quantity alone. The $ T is a very useful input 
quantity for the NN because the W — > Iv decay in the 
signal yields a $ T distribution with higher values than 
the backgrounds. Similarly, the azimuthal angle distri- 
bution between the W lepton and the $J T is useful for 
distinguishing WZ from the backgrounds because they 
do not contain W decays. The total energy transverse to 
the beamline deposited by the WZ decay compared to 
that of background processes and lepton flavor combina- 
tions (eee, ee\x, e/i track, etc.) are also examples of NN 
input variables used. Figure Q] shows the output of the 
NN treatment, with background-like events in simulation 
and data trending toward a value of —1 while signal- like 
events trend towards +1. Note that tt is represented in 
Fig. [TJ but has too small of a contribution to be visible. 

The measured cross section for W Z is extracted from 
the NN output in Fig. [T]with a binned maximum likeli- 
hood fit method. The likelihood function is formed from 
a product of Poisson probabilities for each bin in the NN 
output and Gaussian constraints are applied correspond- 
ing to each systematic uncertainty: 



C 



n 



(i) 



where /ii is the total expectation in the z-th bin, is the 
number of data events in the z-th bin, and S c is a floating 
parameter associated with the systematic uncertainty c. 
The fii is given by 



6 



WZ Signal (Z-Peak Region): All Jet Bins 
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FIG. 1: The NN output for discriminating the WZ signal 
events from background processes within the selected signal 
sample. Note that the ti contribution is small enough to not 
be visible. 
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k 



.measured 



T-cxpcctcd 



ri( i + ^ 



^ jy-expcctcd^ 



(2) 

The sum k is over processes contributing to events in 
bin i and ft is the fractional uncertainty for the pro- 
cess k due to the systematic uncertainty c. Some sys- 
tematic uncertainties are common to more than one pro- 
cess and so are correlated. These correlations are ac- 
counted for in the definition of /Zj through the /| param- 
eters. The (7\r^ x P cctod ^ j s t ne expected number of events 
from process k in the i-th bin. All the background pro- 
cesses are constrained to their SM expectations by set- 
ting the proportion of measured to expected cross sec- 
tion to unity. The likelihood is then maximized with 
respect to the floating systematic (S c ) and cross section 

proportion ( a e XD cctcd 1 parameters, where ,7 ex P ectcd is 
\ a Jwz 

the expected signal cross section and cr measured [ s the WZ 
cross section ultimately measured from the data. This 
method gives a measured value for the WZ cross section 
of a(pp ->• WZ) = 3.93t^(stat)tjj;^(syst) pb, which 
is in good agreement with the aforementioned standard 
model prediction of 3.50 ± O.210, Q. 

The shape and normalization of the px spectrum of 
the Z boson (Fig. [5]) are used to place limits on anoma- 
lous TGCs. The most general modification of the WW Z 
vertex preserving C and P separately is parametrized by 
X z , 9i and n z I2|. In the SM, X z = Agf = Ak z = 
where Agf and Ak z are used to denote the deviations of 
gf and k z from their SM values. In general, the param- 
eters A z , Agi and Ak z can be functions of the invariant 
mass y/1 of the WZ system. Non-zero values of X z , gf 
and k z at large s violate unitarity. To avoid this, each 
coupling is modified by a form factor a(s) = 



TABLE II: Expected and observed limits on anomalous 
TGCs. For each coupling limit set, the two other couplings 
are constrained at their SM values. Values of the couplings 
outside the given observed range are excluded at the 95% 
confidence level. 





A (TcV) 








Ak z 


Expected 


1.5 


(- 


-0.11, 0.12) (- 


-0.12, 0.23) (- 


-0.58, 0.94) 


Observed 


1.5 


(- 


-0.09, 0.11) (- 


-0.09, 0.22) (- 


-0.42, 0.99) 


Expected 


2.0 


(- 


-0.10, 0.10) (- 


-0.11, 0.20) (- 


-0.53, 0.86) 


Observed 


2.0 


(- 


-0.08, 0.10) (- 


-0.08, 0.20) (- 


-0.39, 0.90) 



where evo is the unmodified coupling X z , gf or k z . 

The likelihood of the Z px distribution for various 
anomalous TGC models is used to set limits. The 
expected Z px distribution for a given TGC before 
the effect of the detector response is obtained using 
MCFM The detector acceptance and efficiency are 
modeled by multiplying the MCFM distribution by a Z 
PT-dependent factor. This factor is calculated using six 
different simulated event samples generated at different 
TGC values with the full detector response simulated by 
GEANT3. The TGC values are chosen to be in the param- 
eter space near the existing limits. For each sample, the 
product of acceptance and efficiency is extracted from 
the simulation as a ratio of the reconstructed and gen- 
erated yields. These ratios are averaged together as a 
function of Z px using the maximum variation as an es- 
timate of the uncertainty due to assuming the efficiency 
and acceptance are not dependent on the TGC values. 

A likelihood for each of the couplings, £(X Z ), £(A<?f ), 
and C(Ak z ), is computed as a product of the Poisson 
probability of each of the bins of the Z px distribution for 
the assumed anomalous coupling. Then 95% confidence 
levels are set where (— 21n£) — (— 2\n£ min ) = (1.96) 2 . 
The systematic uncertainties include everything consid- 
ered for the WZ cross section and the additional px- 
dependent uncertainty on the efficiency described pre- 
viously which ranges from 5% to 20%. Systematic un- 
certainties are implemented in a way which most reduces 
the TGC limit sensitivity when fluctuating the signal and 
background by one standard deviation, thereby taking 
a conservative approach in assigning systematic uncer- 
tainty. The observed 95% confidence level limits are con- 
sistent with expectations as shown in Table [TTJ 

To summarize, the WZ production cross section has 
been measured in pp collisions at ^fs = 1.96 TeV 
from reconstructed events in the trilepton plus $ T fi- 
nal state using a likelihood ratio formed from a Neu- 
roBayes neural network distribution that discriminates 
signal from background. This result, a(pp — > WZ) = 



an 



3.93_ ; 53 (stat)_ Q ; 46 (syst) pb, is the most precise mea- 
surement at this energy with an overall uncertainty of less 
than 20% and in agreement with SM predictions. The 
same event sample is also used to perform the most sen- 
sitive probe to date at this energy of anomalous WW Z 
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FIG. 2: The Z pr distributions for data compared to the 
SM expectation for signal (WZ) and background. Also pre- 
sented is how the signal expectation would change with the 
introduction of anomalous couplings near the observed limits. 



couplings. The Z px distribution of the sample is found 
to be in agreement with the SM expectation and is used 
to place limits on anomalous triple gauge couplings. 
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